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We consider a particle moving on the half line x > and subject to a constant 
force in the —x direction plus a delta-correlated random force. At x = the particle 
is reflected inelastically. The velocities just after and before reflection satisfy vj = 
—rvi, where r is the coefficient of restitution. This simple model is of interest in 
connection with studies of driven granular matter in a gravitational held. With an 
exact analytical approach and simulations we study the steady state distribution 
function P{x, v). 
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I. INTRODUCTION 

Consider a particle moving on the half line x > and subject to a constant force in the 
—X direction, such as gravity, and an additional delta-correlated random force with zero 
mean. The Newtonian equation of motion has the form 

^ = -9 + v{t), {vmt')) = 2AS{t-t'), (1) 

where g is a positive constant. With no loss of generality we set = 1, A = 1 throughout 
this paper, since this can be achieved by rescaling |l| of x and t. At the boundary x = the 
particle is reflected inelastically. The velocities just after and before reflection satisfy 

Vf = -rvi, (2) 

where r is the coefficient of restitution. The simple model defined by Eqs. and Q is of 



3|. Cornell, Swift, and Bray 



interest in connection with driven granular gases, see, e.g., 
originally proposed the model and studied the case g = 0. 

According to Eq. ((T)) the mean energy E{t) = (t)^ + a;(t)) of the randomly accelerated 
particle at time t increases as E{t) = E{0) + 1 between boundary collisions but decreases, 
for r < 1, in the collisions. Eventually an equilibrium is reached. In this paper we study the 
equilibrium distribution function P{x, v) with an exact analytical approach and simulations. 
In related earlier work, but without gravity, we have considered (i) the equilibrium of a 
randomly accelerated particle moving on the finite line between two inelastically refiecting 
boundaries (^51, and (ii) non-equilibrium properties of a randomly accelerated particle on 
the half line^, 3] and on the finite line for absorbing, partially absorbing, and inelastic 
boundary conditions . 

The equilibrium distribution function P{x, v) of a particle, which moves according to Eq. 
(0) for X > and is refiected inelastically at x = 0, satisfies the steady-state Fokker-Planck 
equation j^l 

' d d (9^ \ 

-^-l-^:^]Pi-^-) = ^^ (3) 



dx dv 9f " ^ 
with the boundary condition 

P{0,-v) =r^P{0,rv) , v>0. (4) 
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The factors of r in Eq. Q take the change in velocity © into account and ensure that the 
incident and reflected probability currents at the boundary have equal magnitude 

/•OO /'OO 

/= / dvvP{0,-v)= / dvvP{0,v). (5) 
Jo Jo 



'0 JO 

Here we have assumed the standard normalization 



dx / dvP{x,v) = 1 , (6) 

J-oo 

of the distribution function. 

The main goal of the paper is to solve the steady-state Fokker-Planck equation ^ with 
boundary condition (@]). In Section II we derive an integral equation for the distribution 
function P(0, v) at the boundary and show how to obtain P{x, v) from P(0, v) by integration. 
In Section III the asymptotic form of P(0, v) for large and small v is determined, for arbitrary 
r, from the integral equation. For two particular values, r = ^ and |, the integral equation is 
solved exactly and the corresponding P{x, v) obtained by integration. We have also carried 
out numerical simulations of the randomly accelerated accelerated particle in a gravitational 
field. The algorithm is described in Section IV, and the results of the simulations are 
compared with our analytic predictions. Section V contains closing remarks. 

II. INTEGRAL EQUATION FOR P{x, v) 

The Fokker-Planck equation is readily solved by separation of variables. The most 
general solution which vanishes for x ^ oo and v ±oo is given by ^ 

POO 

P{x,v) = e-^l'' / dFw;(F)e-^^'V'i/4,F(-t;) , (7) 
Jo 

where is an arbitrary function, and iIjs,f{v), with F > 0, is the solution of the ordinary 
differential equation 

+ Fv + s) ^sA^) = (8) 



dv"^ / 
which vanishes for v — > ±oo. This solution is given by 

^^^p(v) = F-'/^Ai {F'^h + F-^/h) , (9) 

where A\{z) is the standard Airy function lio|. Normalized as in Eq. the V's.fI^) have 
the orthogonality and closure properties j3| 

, upper sign 

dvvi)s,F{^^)i^s,G{-'") = \ (10) 

6{F — G) , lower sign 
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dF [ i)s,F{-v)^ps,F{-u) - i)s,F{v)i)s,F{u) ] = V ^6{v - u) . (11) 

Makina; use of Eqs. dZj), ©, and the relation rft e-''*Ai(-t) = exp (a^/3) for a > 0, 
one finds that the probabihty density P{x) for the position of the particle is given by 

/OO /"OO 
dvP{x,v)= / dF F-^^^w{F)exp[-{12F)-'^ - Fx] , (12) 
OO Jo 

and that the normalization condition © is satisfied if 

/»oo /»oo 

/ dxP{x)= dFF-^^^w{F)exp[-{12F)-^] =1 . (13) 
io Jo 

To invert Eq. (jlj), it is convenient to use the auxiliary functions ^ 

<PsA^) = ^sA^) - ^ — i'sM-v) > (14) 

which vanish for f > and have the biorthogonality property 

dvv(PsA-^)M-^) = ^{F - G) (15) 



Jo 



on the semi-infinite interval v > 0. Functions of this type were first constructed by Marshall 
a„d Wat.o„ a a„d exp.o.ted in fly Q in ca.cu.a«ng fl.. passage p..ope.*.. 
Using property (fT5|l to solve Eq. for w{F), we obtain 

poo 

w{F)= dvv<pi/iA-vV'^P{^^v) ■ (16) 
Jo 

Substituting this result back in Eq. ((Zj) and imposing the boundary condition (jH) leads to 

/■OO 

r^P{0,rv)= duuGg{v,u)P{0,u) , (17) 
Jo 

POO 

Ggiv,u) = e^^+^^Z^ / dF^,/,A^)(l),/,A-u) ■ (18) 



Equations ^ and (fTBjl - (fTHj) provide us with a strategy for calculating P{x, v). First P(0, v) 
is determined by solving the integral equation (fTTj). Then w{F) and P{x,v) are obtained 
by evaluating the integrals in Eqs. (fTBj) and ((Tj), respectively. 

The kernel in integral equation ()17|1 has a simple physical interpretation. The quantity 
vGg{v, u)dv is the probability that a particle moving according to Eq. with initial position 
a; = and initial velocity u > arrives with speed between v and v + dv on its first return 
to X = 0. Integral equation ()17|) for the equilibrium distribution follows immediately from 
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this interpretation. For consistency with the probabihstic interpretation, Gg{v,u) should be 
symmetric in v and u and normahzed so that 

dvvGg{v,u) = l . (19) 







That the expression for Gg{v,u) in Eq. (fTHj) satisfies both these conditions can be shown 
with the help of Eqs. (HH), and JTH). 

The integral representation of Gg{v,u) in Eqs. (fTHjl is useful because of the simple or- 
thogonality and closure properties of the basis functions. Two other useful representations, 

Jv+u)/2 roo sinhfvr^) 
G,{v,u) = / ^77 \/ ]^^^ K,,{v)K,,{u) (20) 

n^vu Jq cosh(7r7/3) 

Gg{v, u) = A(^^)i/2e{-+«)/2 r I gxp [-(t;2 _ yu + + 3i;V)'^'] 

271" Jq 

X [{v^ - vu + u"^ + 3vuy^)-^^^ + {v^ - vu + u"^ + Svuy^y^]} , (21) 



the standard modified Bessel function 
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follow from a classic paper of McKean 13|, as shown in Appendix A. In Eq. (pUj) . K,^{z) is 



14| . Expressions (plj) and (PT|) are particularly 



convenient for analytical (see Appendix B) and numerical calculations, respectively. 

III. SOLUTION OF INTEGRAL EQUATION (|T7jl 

A. Asymptotic form of P{0, v) for small v 

First we derive the asymptotic form of P(0, v) for small and large v from integral equation 
fjl7|) . Using representation ()21|) for the kernel, one finds that the ansatz 

P(0,t;) , v\0 (22) 

is consistent with the integral equation if the exponent (3{r) satisfies 



/ sm I TT 



l/{/3-2) 

(23) 



The dependence of f3 on the coefficient of restitution r is shown in Fig. 1. As r decreases 
from 1 to 0, (3{r) increases monotonically from to 5/2. The probability density for finding 
the particle at the boundary with velocity v diverges as v approaches zero, and the greater 
the inelasticity, the stronger the divergence. 
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The exponent (3 has the value 2 at 

'Vv^ = 0.163. . . (24) 



e 



For r < Tc, (3{r) > 2, so that P{0,v) diverges more strongly than in the limit f \ 0. 
Thus, for r < Tc the mean collision rate, defined by the integral (0), diverges at the lower 
limit, i.e., the randomly accelerated inelastic particle makes an infinite number of collisions 
in a finite time. Cornell et al. 0] argued that even in the absence of a gravitational field 
this is the case. 

The small v behavior ()22|1 , ()23|) is precisely the same as for a randomly accelerated particle 
moving on the finite line < a; < 1 with inelastic collisions at a; = and 1 and with no 
g.av,tat,o„al field QQ. I„ both case. tKe d.ve.ge.t behavo. fo. . ^ co,„e. .on, the high 
rate of low velocity inelastic collisions with one of the walls. For short times the diffusive 
change in velocity due to the random force, which grows as t^^^, exceeds the change due 
to the gravitational force, which is proportional t. Thus, the random force is primarily 
responsible for the rapid return to the boundary and determines the small v behavior. 

In both cases - the randomly accelerated particle on the half line with a gravitational 
field considered here, and on the finite line with no gravitational field considered in ^ - the 
distribution function P{x, v) remains smooth, extended, and normalizable for r < Tc. In 
neither case does it collapse on the boundary. Thus, the infinite equilibrium collision rate 
does not lead to localization of the particle at the boundary. 



B. Asymptotic form of P{0,v) for large v 

For large v one finds that the ansatz 

P(0, v) ~ t;— Me-«W" , v^oo. (25) 

is consistent with integral equation ()17j) with kernel ()2H1 if the quantities m(r) and q{r) are 

chosen as follows: For | < r < 1, 

, , ^ Inr 3 1 — r , . 

mir) = 1 — , r / N / 1 , Qv) = ^ , (26) 

^ ^ ln[a;(r)/r] ' ^ 2 1 - r + r2 ' ^ ^ 

where 

1 - (2 
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and for r < ^, 

3 1 

m{r) = - , g(r) = — . (28) 

The functions m(r) and q{r) are shown in Figs. 1 and 2, respectively. As r decreases 
from 1 to |, m(r) increases from to 1 and then changes discontinuously to | for r < |. 
As r decreases from 1 to 0, g(r) increases monotonically from to oo. The functional form 
of g(r) changes at r = i, in accordance with Eqs. ()26|) and (j^Hj) . Both q and dq/dr are 
continuous at r = 1/2, but d'^q/dr^ is discontinuous. 

The functions m(r) and g(r) are both non- analytic at r = i. Since they only characterize 
P{0,v) in the large v limit (see Eq. ((23)); it does not follow that P{x,v) is non-analytic in 
r for fixed finite x and v. We have not found any indication of a phase transition, i.e., a 
non-analyticity of the distribution function P{x,v), as r is varied with x and f fixed. 



r. 



C. Exact P(0,?;) for r = ^ and ^ 

The curves for /3(r) and m(r) in Fig. 1 cross at r = | and |. For these two values of 
the exponents in the asymptotic forms (j^^ and ()25j) for small and large v are the same, 
with /5(|) = m(|) = 1 and /5(|) = m{^) = |. In fact, the functions 

turn out to be exact solutions of integral equation (jl7|) for all < f < cxo, as we discovered 
numerically and prove analytically in Appendix B. The quantities y4i/2 and Ai/a in Eq. 
are normalization constants, which we now determine. 

Substituting Eq. (PIJ|) in Eq. (fTBj) and making use of Eqs. © and (HH), we find that the 
corresponding w{F) are given by 

''"^1/24 17-1/2 ^nol7N-U ». - 1 

3 ' 



(^)^^^A/3F-^/^i^i/6((12F)-i) , r 



where Ki,{z), as in Eq. (pUj). is a modified Bessel function. Together with Eq. (fT^. these 
w{F) imply 

1 / 27 \ 

^"^=3- (31) 

for the normalization constants in Eqs. ()29|) and ((301) • 
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Combining Eqs. (fT^ . (I^Ujl. and (jHTjl . one obtains relatively simple expressions 

. iiro(\/2x/3) , r = i 

P(a;) = 7 3 oiv / ; ' 2 ' 

^/,°°rfFF-iiri/6((12F)-i)exp[-(12F)-i-Fx] , r = |, 

for the probability distribution of position. In the limit x ^ 0, P{x) diverges as logx and 
3;-i/6 fQj, r = ^ and ^, respectively, and as x^^~^^^^ for general r, as in Ref. |5j. In the 



limit X — s> 00, P{x) decays as exp(— A/2a;/3) and exp(— A/2a;/3) for r = ^ and |, 
respectively. 

The distribution function P{x, v) may be calculated numerically for r = ^ and | by 
combining Eqs. ((Tj), ©, (jHUj) . and (jHTj) and carrying out the integration over F numerically. 
The V dependence of P{x, v) for r = ^ and several values of x is shown in Fig. 3. For r = ^ 
the curves are qualitatively similar. In both cases the curve for x = diverges at f = 0, 
in accordance with Eqs. (|22j) and (j23|) . For x > the divergence is replaced by a smooth 
peak, which becomes lower and broader as x increases. 

For large x the height of the peak in Fig. 3 varies as P(x,0) ~ x~^^'^ exp(— a/2x/3) for 
r = i and as P{x,0) ~ x~^/^ exp(— a/2x/3) for r = |. For small x, P(x, 0) ~ x~^/^ and 
x^^^'^ for r = ^ and |, respectively, and P(x, 0) ~ x~^^'^ for general r, as in Ref. 

The curves in Fig. 3 are characterized by the moments 

POO 

(y")^ = p(x)-^ / dv t;"P(x,w) , (33) 



where the normalizing denominator P(x) is defined by Eq. ()12j). Taking the derivative of 
Eq. with respect to x, making use of the Fokker-Planck equation Q, and integrating 
partially with respect to v, one finds that 

{v)cc = , (34) 
{v'), = P{x)-^ / dx P(x) . (35) 



According to Eqs. (0), ()33p. and the net probability current at a distance x from the 
boundary vanishes, as expected in equilibrium. Inserting the asymptotic form of P(x), given 
just below Eq. in Eq. one finds that the root-mean-square width (f^)^^ of the 

peak in Fig. 3 increases as x^/^ in the large x limit, for both r = ^ and |. 
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IV. SIMULATIONS 
A. Algorithm 

In our simulations the motion of the particle is governed by the difference equations 

1 fi . V'"^ 

Xn+l = Xn + WnA„+i - - A^^^^ + f - A^^_^i j + r„+i) , (36) 

Vn+i = Vn- A„+i + (2A„+i)^/^ r„+i . (37) 

Here Xn and f„ are the position and velocity at time t„, and A„+i = tn+i — tn- The 
quantities r„ and s„ are independent Gaussian random numbers with (r„) = = and 
(r^) = {sD = 1. This is the same as the algorithm in Refs. except that a gravitational 

field g = 1 has been included. 

The algorithm ()36|1 . (j37j) is consistent with the exact probability distribution 
Pg'^'^'^{x,v; xo,vo;t) in free space, i.e. in the absence of boundaries, of a randomly accel- 
erated particle in a gravitational field g = 1 with initial conditions xq, vq at t = 0. This 
quantity satisfies the time-dependent Fokker-Planck equation ()A1|) (see Appendix A) with 
initial condition ()A2|) . The exact form of Pj'''"^(x, f ; Xq, fo; ^) follows immediately from Eq. 
flA4|) and the expression for PQ'^'^{x,v;xo,VQ;t) in Eq. (27) of 0] 

In free space there is no time-step error in the algorithm (j36|l . (j37j) . The time step A^+i is 
arbitrary. However, for a particle on the half line a; > 0, trajectories are not generated with 
the correct probability close to the boundary. This is because Pg'^'^^{x,v; xo,vo;t) includes 
trajectories which travel from positive to negative x and return during the time t, violating 
the restriction x > 0. As in Refs. 0,0], we make the time step smaller near the boundary 

n 

to exclude these spurious trajectories. Generalizing the discussion in Ref. 8|] to include 
a gravitational field g = 1, we find that the spurious trajectories occur with negligible 
probability if the time step satisfies 

Xn + - ^Al^, - cAl% > , (38) 

where the constant c is about 5 or larger. 

The most efficient time step is the largest A„+i consistent with inequality (jSHI)j which we 
denote by D(x„,f„). Using a smaller time step slows the simulation without improving the 
accuracy. A convenient choice is A^+i = D(x„,f„) + 6, where a small minimum time step 
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has been included. Without it, the step size decreases to zero as the particle approaches the 
boundary, and it never gets there. 

For small r the typical speed of the particle becomes extremely small after many boundary 
collisions. To simulate the behavior reliably, we reduced the minimum time step 6 as the 

n 

speed decreased, following Ref. JS|- After each boundary collision we set the root-mean 
square velocity change Av = (25)^/^, corresponding to the minimum time step, equal to to 
1/500 of the velocity just after the collision and then used this value of 6 until the next 
boundary collision. This Av is the smallest velocity the algorithm can handle reliably. 



B. Results 

According to Eq. (0) the equilibrium rate of boundary collisions in which the particle is 
reflected with velocity between v and v + dv is given by P(0, v) vdv. We used this relation 
to determine P(0,f) from our simulations. 

In Fig. 4 the simulation results for r = 0.75 and r = 0.25 are compared with the exact 
asymptotic forms (jH)), and for small and large respectively. Varying the 

proportionality constants in Eqs. (j^^ and ()25p shifts the dashed and solid curves vertically 
in the figure, without changing the slope. The proportionality constants were chosen to best 
fit the simulation data, and the fit is excellent. 

In Fig. 5 the simulation results for = | and r = | are compared with the exact 
analytical predictions ()29p and (jHTj) . No adjustable parameters are involved. The agreement 
is excellent. 



V. CONCLUDING REMARKS 



In this paper we have studied the equilibrium of a particle on the half line subject to 
both a random force with zero mean and a constant gravitational force which drives the 
particle toward a single inelastic boundary. It is interesting to compare with the case 
of a randomly accelerated particle moving on the finite line between two inelastic boundaries 
with no gravitational field. In both cases the equilibrium distribution function has the same 
divergent behavior P(0, f ) ~ for small f, given in Eqs. (j^^ . (j^Hj) . As discussed below 
Eq. (jnH), this comes from a high rate of low velocity inelastic collisions with one of the 



11 



boundaries. In both cases the random force, not gravity, is responsible for the rapid return 
of a low velocity particle to the boundary. 

In contrast, the equilibrium behavior for high velocities is different in the two cases. With 
a single boundary and a gravitational field, P(0, v) ~ ■y~™'e~''^ for large v, according to Eqs. 
fj^ - PHjl . With two boundaries and no gravitation, P(0,f) ~ exp[— (t>/t>ch)^], according 
to Eqs. (14) and (16) of Ref. j5|]. Thus, in the case of two boundaries high velocities are 
far more strongly suppressed. This is not surprising. Since a high velocity particle quickly 
travels between the two boundaries with negligible change in velocity, the time between 
inelastic collisions decreases as |f as |f | increases. In the case of a single boundary and 
a gravitational field, on the other hand, a particle which leaves the boundary with a high 
velocity travels far from the boundary before it returns, and the mean time between the 
inelastic collisions increases linearly with |f |. 

We conclude with quantitative answers to some simple questions one might ask on reading 
the title of the paper: What is the mean height (x) of the ball above the inelastically reflecting 
surface, how large is the root-mean-square fluctuation Ax = {{x — (x))^)^/^, and what is the 
equilibrium collision rate J ? In the two cases r = ^ and | where we have found the exact 
equilibrium distribution function P{x,v), Eqs. (0), (j^ . (jSU, and (jH^ imply 

r = i : (x) = 6 , Ax = 10.4 , / = i , 

r = 1 : (x) = 3.89 , Ax = 7.68 , / = f , 

For both values of r the fluctuations in the height of the ball are larger than the mean height. 
For r = |, the mean height is smaller and the collision rate is larger than for r = ^. The 
mean height is expected to decrease as r decreases but remain positive for all < r < 1, 
since P(x) remains positive and integrable. The collision rate is expected to increase as r 
decreases and is infinite for r < r^, as discussed below Eq. fl24|) . Finally, we note that these 
results for = A = 1 in Eq. ()39|) are easily generalized to arbitrary values of the parameters 

n 

g and A in Eq. (P). In accordance with footnote [1], the numerical values in Eq. (jHIJjl for 
(x) and Ax are multiplied by g~^A^ and for / by g'^Ar^. 
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APPENDIX A: ALTERNATE EXPRESSIONS FOR Gg{v,u) 

To discuss first-passage properties, we define Pg{x, v;xo,Vo; t)dxdv as the probability that 
the position and velocity of a particle, moving on the half line x > according to Eq. 
with (7 = A = 1, have evolved from xq, fo to values between x and x + dx, v and v + dv 
in a time t without ever reaching a; = 0. The quantity Pg satisfies the time-dependent 
Fokker-Planck equation 

d d d d"^ \ 
with the initial condition 

Pg{x,v; xo,vo;0) = S{x - xo)S{v - vq) (A2) 

and the absorbing boundary condition 

Pg{0,v;xo,Vo;t) = , v>0. (A3) 

In the absence of the gravitational field, i.e. in the case g = instead of g = 1, the corre- 
sponding probability distribution Po{x, v]Xq,Vq] t) satisfies the same Fokker-Planck equation, 
boundary condition, and initial condition, except that the third term —dPj dv in Eq. (jAljl 
is absent. This leads to the simple relation 

Pg{x, V- xo, v^- 1) = e-(--o)/2-V4 p^^^^ ^. (A4) 

between the distribution functions with and without gravity. 

As discussed below Eq. (fT^ . the quantity vGg{v, u)dv in integral equation (fTTjl represents 
the probability that a particle moving according to Eq. (0) with initial position x = and 
initial velocity u > arrives with speed between v and v + dv on its first return to x = 0. 
Since the return occurs between t and t + dt with probability vPg{0, —v; 0, u; t) dv dt, 

POO 

Gg{v,u)= / dtPg{0,-v;0,u;t) = e^''+^^/^Po{0,-v;0,u;l) , (A5) 
Jo 

where we have made use of Eq. ()A4|) and introduced the Laplace transform 

POO 

Po{x,v;xo,Vo;s) = dte~''^Po{x,v;xo,Vo;t) . (A6) 







assic paper on the first passage properties of a randomly accelerated particle, 
I3I derived the exact form of Po(0, — "w; 0, m; t) and its Laplace transform. The 



In a c 
McKean 

expressions for Gg{v,u) in Eqs. and (PT|) follow from Eq. ()A5|) and McKean's Eqs. (5) 
and (6). 
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APPENDIX B: DERIVATION OF EXACT P(0, v) FOR r = | AND i 

Expression for Gg{v,u) and the Kontorovich-Lebedev transforms [3| 

/•OO 

Jo 

POO 

^(7) = 271-^7 sinh(7r7) / dvv~^ f{v)K^^{v), (B2) 

imply 

w-^e^/2^i^(^) = 2 cosh (^^) ^ rfMG<,(t;,M)e-"/2;^i^(M) . (B3) 

For Z7 ^ |, Ki^{v) const x t>~-'^/^e~'", and integral equation ()B3|) is entirely equivalent to 
Eq. (fT7|). with = I and P(0,f) oc ^;-3/2g-3i)/2_ 'jj-^^g^ have found the exact P{0,v) for 
r = |, already announced in Eq. ((221) • 

The other result in Eq. ^I^ . that P(0, f ) oc t;~^e~'' for r = | is an exact solution of the 
integral equation ()17|1 . is readily proved by integrating Eq. ()B3j) over 7 from 7 = to cxo 
and then using the results, 

" cosh (^) i^,,(tO = I _ 

These relations, which are examples of the Kontorovich-Lebedev transforms (jBljl . ()B2jl . 
follow directly from Ref. Q, Eq. (6.795-1). 
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Figure Captions 

Figure 1: The exponents m(r) and /3(r), given by Eqs. ((221), (ESI), and (I23)-(I2HI)- Note the 
discontinuity in m(r) at r = |. 

Figure 2: The decay constant g(r), given by Eqs. (pUj) and (|^. Both g(r) and dq{r)/dr 
are continuous at r = |, but cPq/dr'^ is discontinuous. 

Figure 3: Exact P{x,v), given by Eqs. ((Tj), (fHUjl . and (|^. for r = ^ and, from top to 
bottom, a; = 0, 0.001, 0.01, and 0.1. 

Figure 4: Double-logarithmic plot (base 10) of P{0,v) for r = 0.75 and 0.25. The points 
are the results of our simulations. The dashed and solid lines show the predicted 
asymptotic forms (j^ . for small v and (j^ - lEH)) for large v, respectively, with 
proportionality constants chosen to fit the simulation data. 

Figure 5: Double-logarithmic plot (base 10) of P{0,v) for r = ^ and |. The points are 
the results of our simulations. The solid curves show the exact theoretical predictions 
and (jSH), with no adjustable parameters. 
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FIG. 1: The exponents m(r) and /3(r), given by Eqs. (dH), and (dSJ-llHI)- Note the disconti- 
nuity in m{r) at r = ^. 
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FIG. 2: The decay constant q{r), given by Eqs. H26() and (P5|) . Both g(r) and dq{r)/dr are 
continuous at r = i, but Sqjdr^ is discontinuous. 
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FIG. 3: Exact P{x,v), given by Eqs. 0, ®, (EOI, and (jSB) for r = ^ and, from top to bottom, 
X = 0, 0.001, 0.01, and 0.1. 
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FIG. 4: Double-logarithmic plot (base 10) of P{0, v) for r = 0.75 and 0.25. The points are the 
results of our simulations. The dashed and solid lines show the predicted asymptotic forms (|22|) . 
(^5]) for small v and H25() - ((^ for large v, respectively, with proportionality constants chosen to fit 
the simulation data. 
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FIG. 5: Double-logarithmic plot (base 10) of P{0,v) for r = ^ and 3. The points are the results 
of our simulations. The solid curves show the exact theoretical predictions H29|) and with no 
adjustable parameters. 



